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Abstract: Recently, interesting braneworld cosmologies in the Randall-Sundrum scenario 
have been constructed using a bulk spacetime which corresponds to a charged AdS black 
hole. In particular, these solutions appear to 'bounce', making a smooth transition from a 
contracting to an expanding phase. By considering the spacetime geometry more carefully, 
we demonstrate that generically in these solutions the brane will encounter a singularity in 
the transition region. 
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1. Introduction 

Braneworlds [1, 2, 3] have generated enormous interest in higher-dimensional spacetimes 
amongst particle theorists. A key ingredient in these brane models is that the Standard 
Model particles remain confined to a (3+l)-dimensonal brane, while only the gravitational 
excitations propagate through the full spacetime. Such scenarios provide a new framework in 
which to address many longstanding puzzles in particle physics, such as the hierarchy problem. 
The cosmology community has also shown an increasing interest in braneworlds [4, 5], since 
this is another field where brane models have the potential to provide novel solutions to many 
of the perennial questions. 

In the present paper, we will focus on one small aspect of the braneworld description 
of cosmology. In particular, we are interested in a certain family of cosmological solutions 
[6] which were recently proposed in the context of the Randall-Sundrum (RS) scenario [2, 
3]. Recall that the RS model introduces a codimension-one brane into a five-dimensional 
bulk spacetime with a negative cosmological constant. The gravitational back-reaction due 
to the brane results in gravitational warping which produces massless graviton excitations 
localized near the brane. Fine tuning of the brane tension allows the effective four-dimensional 
cosmological constant to be zero (or nearly zero). Brane cosmologies where the evolution is 
essentially that of a four-dimensional Friedmann-Robertson- Walker (FRW) universe can be 
constructed with a brane embedded in cither AdS [4] or an AdS black hole [5, 7]. 

In either of the above cases, however, the cosmological evolution on the brane is mod- 
ified at small scales. In particular, if the bulk space is taken to be an AdS black hole with 
charge, the universe can 'bounce' [6]. That is, the brane makes a smooth transition from a 
contracting phase to an expanding phase. From a four-dimensional point of view, singularity 
theorems [8] suggest that such a bounce cannot occur as long as certain energy conditions 
apply. Hence, a key ingredient in producing the bounce is the fact that the bulk geometry 
may contribute a negative energy density to the effective stress-energy on the brane [9]. At 
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first sight these bouncing braneworlds are quite remarkable, since they provide a context in 
which the evolution evades any cosmological singularities yet the dynamics are still controlled 
by a simple (orthodox) effective action. In particular, it seems that one can perform reli- 
able calculations without deliberating on the effects of quantum gravity or the details of the 
ultimate underlying theory. Hence, several authors [10, 11, 12] have pursued further develop- 
ments for these bouncing braneworlds. In particular, ref. [12] presents a critical examination 
of the phenomenology of these cosmologies. 

In the following we re-examine these bouncing brane cosmologies, paying careful attention 
to the global structure of the bulk spacetime. We find that generically these cosmologies are 
in fact singular. In particular, we show that a bouncing brane must cross the Cauchy horizon 
in the bulk space. However, the latter surface is unstable when arbitrarily small excitations 
are introduced in the bulk spacetime. The remainder of the paper is organized as follows: 
We review the construction of the bouncing braneworld cosmologies in section 2. Section 
3 presents a discussion of the global structure of the full five-dimensional spacetime and 
the instability associated with the Cauchy horizon. We conclude in section 4 with a brief 
discussion of our results. 



2. Construction of a Bouncing Braneworld 

We consider a four-dimensional brane coupled to five-dimensional gravity with the following 
action 

'-TdalL '^^^ h + i - - si^ L ^^^^ + //'^^^ ■ 

Here, i?5 denotes the Ricci scalar for the bulk metric, g^^^,, and F^i, is the field strength of a 
bulk gauge field. The (negative) bulk cosmological constant is given by As = while 
the brane tension is T = ^^q^^ ■ The length scales L and A are introduced here to simplify the 
following analysis. The induced metric on the brane is denoted by jab- With the last term in 
the action (2.1), we have allowed for the contribution of extra field degrees of freedom which 
are confined to the brane, e.g., the Standard Model fields in a RS2 scenario [3]. 

The bulk equations of motion are satisfied by the five-dimensional charged AdS black 
hole solution with metric 

dsl = -V{r)dt' + y^^+r'di:l, (2.2) 

where 

no-^+fc-^ + l^ (2.3) 

and the gauge potential is At = In the metric above, dS| denotes the line element on 
a three-dimensional sphere, flat space or hyperbolic plane for k = +1, or —1, respectively 
(with unit curvature for the cases k = ±1). The parameters n and q appearing in the solution 
are related to the ADM mass and charge of the black hole — see, e.g., [12, 13]. Note that 
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this solution contains a curvature singularity at r = 0, but if fx is large enough, there are 
two horizons at radii r = r± solving V{r±) = 0. A Penrose diagram illustrating the maximal 
analytic extension of such a black hole spacetime is given in Figure 1. In different parameter 
regimes, the positions of these two horizons may coincide (or vanish, i.e., r± become complex) 
to produce an extremal black hole (or a naked singularity). We will not consider these cases 
in the following. 

The brane is modelled in the usual thin-brane 
approximation. That is, its world volume is a hy- 
persurface, B, which divides the bulk spacetime, 
Ai, into two regions. At this hypersurface, the 
bulk metric is continuous but not differentiable. 
Using the standard Israel junction conditions [14] 
(see also [15]), the discontinuity in the extrinsic r = 
curvature is interpreted (5-function source of 
stress-energy due to the brane. Then, defining the 
discontinuity in the extrinsic curvature across B as 



/Cafe = K, 

by 



ab 



K the surface stress-tensor is given 



Sab 



1 



{K^ab - lab ^"c) 



(2.4) 




Figure 1: Penrose diagram for maximally 
extended AdS Reissner-Nordstrom black 
hole. 



In the case of an empty brane with only tension 
(i.e., a brane on which no internal degrees of free- 
dom are excited), one has Sab = —T^ab- 

The construction of the braneworld cosmology 
[7] then proceeds by taking two copies of the AdS 
black hole geometry, identifying a four-dimensional 
hypersurface r = a(r), t = h{T) in each, cutting out 
the spacteime regions beyond these hypersurfaces 
and 'gluing' the two remaining spacetimes along 
these surfaces. While asymmetric constructions are possible (see e.g., [12]), we will focus on 
the case where the two bulk spacetime geometries are characterized by the same physical 
parameters {fi,q,L). With this choice, the calculation of the surface stress-tensor simplifies, 
since ICab = Note, however, that the gauge fields are chosen with opposite signs on 

either side of the brane. Then the flux lines of the bulk gauge field extend continuously 
over the brane, starting from a positively-charged black hole on one side and ending on the 
negatively-charged one on the other. In this case, the brane carries no additional charges. 
We will return to consider a charged brane in the discussion section. Since the black hole 
geometry includes two separate, asymptotically AdS regions, an economical approach to this 
construction would be to glue together two mirror surfaces in each of the asymptotic regions.^ 

Of course, the hypersurface described above must be determined to consistently solve 
the Einstein equations (or alternative lv. the Israel junction conditions (2.4)) for a physically 
^Note that this periodic construction is distinct from the RSI models [2], e.g., there is a single positive 
tension brane here, rather than two branes one of which has a negative tension. 
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reasonable surface stress-tensor. Here we follow the analysis of ref. [7]. Identifying the time 
coordinate on the brane as the proper time, r, fixes 



F(a)6^ = j|^ + l. (2.5) 

The induced metric then takes a standard FRW form: 

ds^ = lah dx^dx^ = -dr'^ + a(r)^dS^ . (2.6) 

Again, the brane worldvolume in the bulk spacetime (2.2) is given by r = a(r) and t = 5(r) 
and so the Israel junction conditions (2.4) imply 

4^^, (2.7) 

a A3 

where the 'dot' denotes dr, and we have included a homogenous energy density p for brane 
matter. Stress-energy conservation would imply that the latter satisfies p + 3^ {p + p) = 0, 
where p is the pressure due to brane matter. 

A conventional cosmological or FRW constraint equation for the evolution of the brane 
is produced by squaring eq. (2.7): 

Here, we have defined a 'vacuum' curvature scale on the brane as 

Implicitly, i'^ is assumed to be positive here, which leads to the cosmological evolution being 
asymptotically de Sitter. However, this assumption is inconsequential for analysis of the 
cosmological bounce which follows. We can also write out the effective cosmological and 
Newton constants in the four-dimensional braneworld as 

where the latter comes from matching the term in eq. (2.8) linear in p to the conventional 
FRW equation in four dimensions: (|) + ^ = ^p. Of course, the FRW constraint in this 
braneworld context also comes with an unconventional term quadratic in p [4]. 

The bulk geometry introduces various sources important in the cosmological evolution 
of the brane. The mass term, /x/a^, behaves like a conventional contribution coming from 
massless radiation. The charge term, —q^/a^, introduces a more exotic contribution with 
a negative energy density. This is another example of the often-noted result that the bulk 
contributions to the effective stress-energy on the brane [16] may be negative — see, e.g., [9]. 
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Many exact and numerical solutions for the Friedmann equation (2.8) can be obtained in 
various situations, e.g., [10, 11]. However, one gains a qualitative intuition for the solutions 
in general by rewriting eq. (2.8) in the following form: 



= + U{a) , 



(2.11) 



where 



U{a) = V{a) - 



A2 



(2.12) 



a2 ^ 3a4 £2 



U(a) =V(a) - |, 



and for simplicity we have assumed an empty brane, i.e., p = 0. In this form, we recognize 
the evolution equation as the Hamiltonian constraint for a classical particle with zero energy, 
moving in an effective potential U (a). In this case, the transition regions where the braneworld 
cosmology 'bounces' are identified with the turning points of the effective potential. We have 
also expressed the latter in terms of the metric function V{a) in eq. (2.12) because we will 
want to discuss the position of the turning points relative to the position of the 'horizons', i.e., 
r± . Recall that we assume the bulk solution corresponds to a black hole with a nondegenerate 
event horizon. That is, wc will assume that there arc two distinct solutions, r±, to V{r) = 0. 
Then, there are two physically distinct possibilities for a bounce. 

The first only occurs with k = +1, i.e., 
with a spherical brane world, and positive £^ 
(or equivalently A4 > 0). In this case, at 
large a, the effective potential becomes large 
and negative. The next most important con- 
tribution at large a is the constant term k 
and hence if A; = +1, the potential may have 
a zero at large a. This bounce is typical of 
those one might find in a de Sitter-like space- 
time, e.g., [17]. It is driven by the spatial 
curvature and occurs as long as the effective 
energy density from the bulk black hole con- 
tributions or braneworld degrees of freedom 
is not too large. The turning point occurs 
at some large ads and in particular, it is not 

difficult to show that a^s > r^. That is, the brane bounces before reaching the black hole. 
In fact the presence of the black hole with or without charge is really irrelevant to this kind 
of bounce. For example, setting /j, = = q in eq. (2.12) produces a de Sitter cosmology on 
the empty brane. 




Figure 2: Effective potential, U{a), appearing in 
eq. (2.11) for the evolution of the scale factor, a(T). 
The turning point, Ob, occurs inside the Cauchy 
horizon r_. 
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The second type of bounce is generic for a wide range of parameters. It occurs at small a 
where the positive ja^ term dominates the potential (2.12), i.e., where the exotic negative 
energy dominates the Friedmann constraint (2.8). As is clear from the first line of eq. (2.12), 
[/(a) < Via) and therefore the turning point occurs at aj,, inside the position of the Cauchy 
horizon, i.e., < r_, as illustrated in Figure 2. The latter result will be essential in the 
following discussion. 

The Penrose diagram for the bouncing braneworld cosmologies is shown in Figure 3. In 
the 'cut and paste' procedure outlined above, the singularity on the right side of the first 
black hole is cut out, but the singularity on the left remains. Hence, the remaining portion 
of the r = r_ surface is still a Cauchy horizon. 



Figure 3: Penrose diagram for a bouncing braneworld model with ingoing modes at the event and 
Cauchy horizons. The grey areas are those regions of spacetime that are cut out in the construction, 
with identification performed along the boundary. 

Note that in Figure 3, the brane trajectory enters the region between the horizons across 
the segment AB and exits across the opposite segment CD. One can verify that this occurs 
in all cases using eqs. (2.7) and (2.5). From the latter, we find that 



If the brane tension and the energy density, p, are both assumed to be always positive, then 
the last factor is always positive. Furthermore, for r_ < a < r_(_, V{a) < 0. Hence, the 
right hand side above is non-zero and has a definite sign for the entire range r_ < a < r-|_. 
Therefore t cannot change sign along the brane trajectory within the black hole interior. It 
then follows that if a trajectory starts at a point on AB with {t,r) = (oo,r^), then it must 



C 





(2.13) 
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run across the black hole interior to a point on CT> with {t,r) = (—00, r_) — see Figures 1 
and 3. 



3. Instability Analysis 

In the previous section, we reviewed the construction of a broad family of bouncing braneworld 
cosmologies [6]. A key result was that the turning point for the brane's trajectory in the bulk 
geometry was inside the Caiichy horizon of the charged AdS black hole. However, previous 
studies in classical general relativity found that the Cauchy horizon is unstable when generic 
perturbations arc introduced for charged black holes in asymptotically flat [18], or de Sitter 
[19, 20] spaces. Below, we will show that the same instability arises in the asymptotically 
AdS case as well. This is problematic for the bouncing braneworld cosmologies, as generically 
the contracting brane will reach a curvature singularity before it begins re-expanding. 

In the following, we demonstrate the instability of the Cauchy horizon to linearized 
perturbations in the bulk. Our approach will be two-fold. We begin by examining linearized 
fluctuations of a massive Klein-Gordon field propagating in the background. Secondly, we 
consider gravitational and electromagnetic perturbations. In both cases, it is found that 
an observer crossing the Cauchy horizon would measure an infinite flux from these modes. 
The expectation is then that the full nonlinear evolution, including the back-reaction on the 
background metric, will produce a curvature singularity. 

Many of the expressions appearing in the linearized analysis involve the surface gravities 
of the two horizons in the background. The surface gravities are given, as usual, by 

dV 



1 



dr 



(3.1) 



=r± 



An important observation in the following is that k^/k^ < 1, which follows from r_ < 
Now it will be convenient to define the event and Cauchy horizons implicitly by re-expressing 
the metric function (2.3) as 

y^^^ ^ (r^-4)(r^-r^)(r^ + rg) 

This expression also defines ro as determining the complex roots of V{r). Further, the analysis 
is facilitated by introducing some new coordinates to describe the background geometry (2.2). 
In particular, we define the radial tortoise coordinate 

j.^ = .J_ log 1^ "^^+1 3_ log l'^"^-! _! ^0-^^ tanh"^ — , (3.3) 

* 2k+ r + r+ 2«;_ r + r- (rg + r^)(ro + r^) 

which is chosen to satisfy dr* = dr/V{r). The focus of the following analysis will be the 
behavior of linearized perturbations in the range r_ < r < r+ {i.e., region II in Figure 1). 
In this region, we have r* — ±00 as r ^ r^. Finally, it will be useful to work with null 
coordinates, 

u = r^c — t , V = r^c + t , (3.4) 
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with which the hne element becomes ds^ = V{r)dudv + r^dS^. 

The massive Klein-Gordon equation in the charged black hole background (2.2) may be 
expanded as 

^ -92$ + V{r)dl^ + -^dr {V^ V{r)) dr^ + ^^1^ - M^^ = , (3.5) 



V{r) \/—g f 

where we write for the Laplacian on the three-dimensional space appearing in the 
line element (2.2). The eigenvalue problems for and V| each have known solutions with 
eigenvalues, say, —uP' and — Hence by separation of variables, the Klein-Gordon equation 
is reduced to a single ordinary differential equation for $(r), 

-dlk + Vir) + + ^ + ^> = u;^ i , (3.6) 



where we have introduced the tortoise coordinate (3.3) and rescaled 4> = r^/^$. 

As we approach the Cauchy horizon, the second term on the left hand side of eq. (3.6) 
vanishes, leading to oscillatory solutions exp(ibia;r^,). Now, the flux seen by an observer 
freely falling across the horizon, with five- velocity ^ is proportional to the square of the 
scalar T = U^df^i^. T then includes a contribution proportional to e'*""9„<I> near the Cauchy 
horizon. Since the solutions of eq. (3.6) are oscillatory as r ^ r_, we have that this term, 
and hence T , diverges. Similar divergences appear in the observed energy density for these 
linearized perturbations, and so the expectation is that when back-reaction is included, the 
metric will develop a curvature singularity. 

Next we proceed to a more rigorous analysis of metric and Maxwell field perturbations, 
following the method of Chandrasekhar and Hartle [18]. We are simply establishing the 
existence of unstable modes and so, for simplicity, we fix = and consider an "axial" 
perturbation of one of the flat space coordinates. However, the extension of this analysis to 
general perturbations and backgrounds is straightforward. 

The unperturbed bulk metric (2.2) is 

ds"" = -V{r) dt^ + ^ + ^{dx^ + dy^ + dz^) , (3.7) 
V (r) L'^ 

where V{r) is as given in eq. (2.3) with k = We now focus on a class of perturbations 
where this metric is modified by replacing 

dz dz + qt{t,r)dt + qr{t,r)dr . (3.8) 

Similarly for the Maxwell field, we introduce perturbations: 5F = {ftz{t, r)dt + frz{ti r)dr) A 
dz. The linearization of the bulk Einstein and Maxwell equations about the background 
solution may be reduced to a single Schrodinger-like equation: 
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In this equation, we have defined F = r^^'^ftz and assumed an e^'"^' dependence for all fields. 
To apply the standard results of scattering theory below, it is important to note that the 
effective potential {i.e., the pre-factor in the second term on the left hand side) is bounded, 
negative and integrable throughout the black hole interior. Further, we note that the effective 
potential vanishes as exp (±K-i-r*) for r* — ^ =^oo. The other components of the perturbation 
are related to F by 



fr, = -driF/r'/^) , 



I 

-2 



Qtr = -^i^^F , (3.10) 

where Qtr = dtQr — drqt- Note that the linearized equations only fix the metric perturbations, 
Qr and qt, up to infinitesimal coordinate transformations of z, but Qtr provides a gauge 
invariant combination which is completely determined. 

We introduce a basis of solutions to eq. (3.9) normalized so that near the event horizon, 
i.e., r* — —GO, 

F{u:,u) ^ e'^'-' , 

F(a;,r,) ^ e-^'^'-* , (3.11) 

representing, respectively, ingoing and outgoing modes. The full solution to (3.9) may then 
be written as 

— W {u)F {uo, r*)+ W (a;) F (a;, r*) e'^^ . (3.12) 
-OO 27r L J 

At present, we are only interested in the ingoing modes, whose profile is determined by W {^)- 
The outgoing modes may be similarly dealt with, but extra analysis would be required to show 
they lead to a divergent flux. We will return to this point near the end of the section. 

We are free to choose any reasonable initial profile for the ingoing modes. However, one 
restriction which we impose on the initial frequency distribution W (w) of ingoing modes is 
that an observer falling across the event horizon at AB measures a finite flux. The flux for 
such an observer contains a term T ~ e^'^+"()„F. Hence considering eq. (3.12), we require 
that W (i^) have at least one pole with Dim(a;) < —kj^. 

The initially-ingoing modes are scattered by the potential in region II, leading to both 
ingoing and outgoing modes at the Cauchy horizon. Scattering theory will be used to relate 
the functions introduced in (3.11) to a second set normalized so that, as — >^ oo, 

F (cc;,rOe-^'^* ^ {A{uo)e'''''' + 5(L^)e-*^^*) e"^'^* 

= ^(c^)e*^" + 5(aj)e-*^^ . (3.13) 

Clearly, the dominant contribution to the flux at the Cauchy horizon results from 

/OO _^ 
(LjujW {u)A{u)e''^'' . (3.14) 
-OO 
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In terms of the Schrodinger equation describing the perturbations, it is the modes that 
are "transmitted" across the potential that constitute this potentially-divergent flux. These 
modes skim along just outside the Cauchy horizon heading towards the brane. 

The integral in (3.14) may be computed by closing the contour in the upper-half-plane. 
Using arguments from [18], with virtually no modification, we find that ujA{u)) is analytic in 
the infinite strip between ibi«;+. For simplicity, we'll further assume that W {uj) is analytic 
in the strip [0,iK+] and that it is non-zero for lo = ik^. With these assumptions, the leading 
term in (3.14) is from the residue of the pole at a; = iK+: 



Since k- > k+, this flux always diverges as u — > oo. Relaxing our assumptions on the 
analyticity of W (w) in the upper half plane could lead to additional divergent contributions 
to the flux, but we will not consider those here. 

Note that the branc and boundary conditions at the brane played no role in the scattering 
analysis above. While the brane will affect the complete scattering of modes inside the event 
horizon, the basic source of the instability is the same piling up of infalling modes on the 
Cauchy horizon found in previous examples [18, 19]. Hence we disregard the details of the 
scattering of modes at the brane, just as the original discussion of the instability for the 
Reissner-Nordstrom black hole [18] ignored the presence of a collapsing star forming the 
black hole. 

However, for completeness, let us briefly discuss the boundary conditions which must 
be imposed on the perturbations at the branc. First, the metric perturbations must be 
matched across the brane surface so that no additional contributions are induced in the surface 
stress-energy (2.4). In particular, the axial perturbation (3.8) considered above induces a 
new K-rz component in the extrinsic curvature, and this component must be continuous 
across the brane. Similarly, continuity is imposed on the Maxwell field strength. More 
precisely, to ensure that no electric charges or currents are implicitly induced on the brane, 
we require that all components n^F^yf^ are continuous, where and are the unit normal 
and any tangent vector to the brane. Finally, since we are working with perturbations to 
the field strength directly, and not the gauge potential, we must demand continuity of the 
tangential components, tlFni,t2, to ensure there are no magnetic charges or currents induced, 
i.e., dF = 0. 

We close this section with a discussion of the initially-outgoing modes defined by the 
distribution W {oj) in eq. (3.12). In Figure 3, we will primarily consider modes entering the 
interior region on the left through the lower portion of AB. In this case, to contribute to the 
instability at the Cauchy horizon CT>, these modes must be reflected by the curvature {i.e., 
by the effective potential in eq. (3.9) to become ingoing. This scattering leads to a different 
analytic structure in eq. (3.13) for A[lv), describing the u-dcpcndcnt modes at the Cauchy 
horizon. In the "un-cut" spacetime with no brane in place, this structure is identical to that 
obtained for the contribution of the initially-ingoing modes to Biui). Of course, inserting the 




(3.15) 
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brane in the black hole interior produces a more complicated scattering problem, the details 
of which would depend on the precise brane trajectory. For example, the outgoing flux would 
receive additional contributions from perturbations transmitted across the brane from the 
right hand side of Figure 3, as well as from initially-ingoing modes which are back-scattered 
by the brane. Wc did not attempt a detailed study of these contributions. 

Now , following the standard analysis with no brane in place, we find that the contribution 
to u>A{uj) from the outgoing modes is analytic in the semi-infinite strip (— iAv_). If we 
assume that W (w) is analytic in the strip (0,iK_), then we would find, upon closing the 
contour in the upper-half-plane, that the contribution to the flux is finite. However, it is 
consistent with the requirement that an observer crossing AC measure a finite flux, to allow 
W (uj) to have poles in the range < 3m{uj) < k_. With such a choice, there will be 
divergent contributions to the flux, provided that the residue of A is non-zero at these poles. 
This effect differs from that discussed above in that the leading contribution to the flux 
comes from a pole in the initial frequency distribution rather than the scattering coefficient 
^4(0;). A similar discussion played an important role in demonstrating the instability of the 
Cauchy horizon of de Sitter-Reissner-Nordstrom black holes over the entire range of physical 
parameters [20]. 

4. Discussion 

One of the most interesting features of the brancworld cosmologies presented in ref. [6] is that, 
while they seem to evade any cosmological singularities, their evolution is still determined by 
a simple effective action, albeit in five dimensions. However, our present analysis indicates 
that instabilities arise in the five-dimensional spacetime, and that the brane will generically 
encounter a curvature singularity before bouncing. The two essential observations leading 
to this result were: i) the turning point for the brane cosmology occurs inside the Cauchy 
horizon of the maximally-extended geometry of the charged AdS black hole and ii) a standard 
analysis within classical general relativity shows that the Cauchy horizon is unstable against 
even small excitations of the bulk fields. Note that from these results we cannot conclude 
that the brane does not bounce, but rather due to the appearance of curvature singularities, 
the evolution can not be reliably studied with the original low energy action (2.1). 

Of course, one may ultimately have reached this conclusion since the full bulk spacetime 
still includes a curvatTirc singularity at r = — see Figure 3. However, while the latter 
remains distant from the brane, those at the Cauchy horizon are of more immediate concern 
as they intersect the brane's trajectory. 

In the discussion of metric and gauge field perturbations in section 3, we fixed k = and 
limited ourselves to modes that depended only on t and r to simplify the discussion. One may 
be concerned by the fact that these modes have infinite extent in the three-dimensional fiat 
space and so we present a brief discussion of the full analysis. Generalizing our results to the 
most general perturbation is straightforward but tedious. For an arbitrary linearized pertur- 
bation, the separation of variables would naturally lead to considering Fourier components 
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in the {x, y, z) directions with a factor exp (in • x). Since we require a superposition of these 
modes for many different n to localize the perturbation, we cannot simply rotate in the flat 
space to remove the dependence on one of the spatial coordinates. Thus the general analysis 
necessarily involves an ansatz for the perturbations dependent on all five coordinates, which, 
of course, requires extending the perturbations to additional components of both the met- 
ric and gauge field. Appropriate linear combinations of these perturbations would decouple, 
giving a set of Schrodinger-like equations, similar to that found above. While the potentials 
in each of these equations is different, there are typically simple relationships between them 
implying relations between the solutions — for further discussion of these relations, see [18]. 
Then it is sufficient to solve only one of the equations, and the analysis, and the results, are 
essentially the same as presented above 

Of course, our preliminary analysis with massive Klein-Gordon modes included all of the 
spatial modes, and further applied for all of the possible values of k, specifying the spatial 
curvature on the brane. In all cases, there was an infinite flux of these modes at the Cauchy 
horizon. While further analysis of the full scattering and boundary conditions would be 
required to make this consideration of fluxes rigorous, the end result would be the same. 
Hence we are confldent that the results for the metric and gauge field perturbations with 
k = also carry over for A; = ±1. 

Recall that, as discussed in section 2, an apparently economical approach to constructing 
these bouncing cosmologies would be to cut and paste along two mirror surfaces in each 
of the separate asymptotically AdS regions of the black hole geometry. In such a periodic 
construction the nature of the singular behavior would be slightly different. As discussed 
around eq. (2.13), the brane trajectory is unidirectional in the coordinate time, t. Hence 
in Figure 3, if a brane enters the event horizon to the right of the bifurcation surface A, 
then it must exit through the Cauchy horizon to the left of V. However, the same result 
requires that a brane trajectory entering to the left of A exits to the right of V. Therefore 
in the periodic construction above, the two mirror trajectories must cross at some point in 
the region r_ < r < r+, as illustrated in Figure 4. Hence, the evolution is singular in that 
the fifth dimension collapses to zero size in a finite proper time. One redeeming feature of 
this collapse is that the curvature remains finite, and hence one might imagine that there is 
a simple continuation of the evolution in which this 'big crunch' is matched onto a 'big bang' 
geometry. Similar collapsing geometries have been a subject of great interest in the string 
theory community recently — see, e.g., [21, 22]. Resolving precisely how the spacetime evolves 
beyond such a 'big crunch' is an extremely difficult question and as yet string theory seems to 
have produced no clear answer. In particular, it seems that these geometries are also subject 
to gravitational instabilities [23] not dissimilar to those found here. In the present context, 
the situation is further complicated as the precise matching procedure for the background 
geometry is obscure. Naively, one might be tempted to continue beyond the collapse point 
£ with the doubly shaded region in Figure 4. However, a closer examination shows that 
the brane would have a negative tension in this geometry. The other natural alternative is 
to match the crunch at £ to the big bang emerging from J^, but the gap in the embedding 
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geometry would seem to complicate any attempts to make this continuation precise. In any 
event, it is clear that once again these knotty questions can not be resolved using the low 
energy action (2.1) alone but, rather, one would have to embed this scenario in some larger 
framework, e.g., string theory. 




Figure 4: Penrose diagram for periodic construction of the braneworld cosmology 

Much of our discussion has focussed on the bouncing cosmologies for an empty four- 
dimensional brane, but the analysis and the results are easily extended to other cases. One 
simple generalization would be for higher dimensional cosmologies following, e.g., [11]. The 
instability found here would also appear in the assymetric constructions discussed in ref. [12]. 

A more interesting generalization to consider is adding matter excitations on the brane. 
As long as the energy density is positive, such matter contributions will not affect the result 
that the brane crosses the Cauchy horizon. At first sight, it would also seem that reasonable 
brane matter cannot prevent the bounce. The negative energy contribution arising from the 
bulk charge is proportional to For a perfect fluid (in four dimensions), this would 

require the stiffest equation of state consistent with causality [24], i.e., p = p. For example, 
a coherently rolling massless scalar field would yield p oc Hence it would seem that 

the term —q'^/a^ would dominate the p contribution coming from brane matter and a bounce 
would be inevitable. However, ref. [12] recently pointed out the contribution in the FRW 
constraint (2.8) can prevent a bounce. In fact, with any equation of state p = wp with 
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w > 0, this contribution can dominate the bulk charge contribution. Hence with a sufficiently 
large initial energy density on the brane, a big crunch results on the brane. This crunch 
corresponds to the brane trajectory falling into the bulk singularity at r = 0. It then follows 
that the brane must cross the Cauchy horizon in this case as well, and we expect singularities 
to develop there with generic initial data. 

More broadly, mirage cosmologies [5] are induced by the motion of a brane in a higher 
dimensional spacetime. A general warning which the present analysis holds for these models 
is that Cauchy horizons are quite generally unstable. Hence if a particular solution involves a 
brane traversing such a surface in the bulk spacetime, one should expect that these cosmologies 
will encounter singularities for generic initial data. 

At this point, wc observe that in the literature much of the discussion of these brane 
cosmologies treats the brane as a fixed point of a Z2 orbifold, rather than making a symmetric 
construction as discussed in section 2. As discussed there, one must flip the sign of the gauge 
potential in the background solution on either side of the brane in order that the brane is 
transparent to field lines. In contrast for a Z2 orbifold, the field lines end on the brane. As 
there is no natural coupling of a one-form potential to a thrcc-brane in five dimensions, the 
model must be extended to include charged matter fields on the brane. One comment is that 
as the action (2.1) does not explicitly include these degrees of freedom or their coupling to 
the Maxwell field, we cannot be sure that the analogous construction to that presented in 
section 2 will yield a consistent solution of all of the degrees of freedom. One might also 
worry that the simplest solutions would have additional instabilities associated with having 
a homogeneous charge distribution throughout the brane. 

Bouncing cosmologies have long been of interest [25]. Much of their appeal lies in their 
potential to provide a calculable framework to describe the origins of the universe. Apart 
from the present discussion, braneworlds and higher dimensions have inspired many attempts 
to model a bouncing cosmology, including: pre-big bang cosmology [26] ; cyclic universes [27] 
based on a Lorentzian orbifold model [21]; braneworld cosmologies induced by cyclic motion 
in more than one extra dimensions [28];^ universes with higher form fluxes [30, 31], which 
are related to S-brane solutions [32]; braneworld cosmologies [33] with an extra internal time 
directions [34]. However, as well as the present model, none of these works has yet provided 
a compelling scheme which is free of pathologies or obstructions to prediction. We may take 
solace from the absence of any simple bounce models in that it appears that understanding 
the early universe and, in particular, the big bang singularity demands that we greatly expand 
our understanding of quantum gravity and string theory. 
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